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TECHNICAL MEMORANDUM 1237 


ON THE MOTIONS OF AN OSCILLATING SYSTEM UNDER 
THE INFLUENCE OF FLEP-FLOP CONTROLS* 

By I. Fldgge-Lotz and. K. Klotter 

PART I.- MOTIONS OF AN OSCILLATOR OF ONE DEGREE 
OF FREEDOM; CONTROL WITH REGULATOR-POSITION 
. CONTROL WITHOUT SHIFTINGS OF REVERSALS. 

Abstract: So-called flip-flop controls (also called "on-off— course 

controls") are frequently preferred to continuous controls 
because of their simple construction. Thus they are used 
also for the steering control of airplanes. Such a body 
possesses — even if one thinks , for instance , only of the 
symmetric longitudinal motion — three degrees of freedom 
so that a study of its motions under the influence of an 
intermittent control is at least lengthy. Thus, It Is 
suggested that an investigation of the basic effect of Buch 
a control first be made on a system with one degree of 

freedom^ . Furthermore, we limit ourselves in the present 
report to the Investigation of an "ideal" control where the 
control surface Immediately obeys the command given by the 
"steering control function". Thus the oscillation properties 
of the control surface and the defects in linkage, sensing 
element, and mixing device are, at first, neglected. As 
long as the deviations from the "ideal" control may be 
neglected in practice, also the motion of the control 


*"Uber Bewegungen eines Schwingers unter dem Einfluss von 
Schwarz-Weiss— Steuerungen." Zentrale fhr vissenschaftllches 
Berichtswes9n der Luftfahrtf orschung des Generalluf tzeugmeisters (ZWB), 
Untersuchungen und Mitteilungen Nr. 1326, Berlin, August 1, 19^3* 

^In carrying out these investigations we were effectively 
supported by study councillor K. Scholz (DVL) and H. F. Hodapp (Vlf S) . 
In the beginning Dr. Bader (DVL) also participated in the investi- 
gations. 
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surface takes place at the heat of the motion ccf the 
principal system. The aim of our investigation is to 
obtain a survey of the influence of the system and control 
coefficients on the damping behavior which is to be 
2 

attained . 

Outline : 1. The Flip-Flop Control (System Equation, Steering-Control 

Equation, Equation of motion, Steering-Control Function, 
Systems of Type A and B, Reversal Points of Type a 
and b) 

2. The Finite Equation of Motion, the Equation for the 
Control Function and Its Derivative 

3* Representation of the Course of Motion. _The Generating 
Motion. The Phase Plane. The Plane of the Reversal 
Values — • 

4. Division of the Plane of the Reversal Points into 
"Intercepts" and "Limit Points" ("Starting Points" and 
"End Points", Respectively) and into Reversal Points a 
and Reversal Points b. Respectively 

5 . Curves of Constant Interval Length vt^ (isochrones) in 

the Plane of the Reversal Values and their Particular 
Cases 

6. Energy Consideration and Periodic Solutions 

7 . Comprehensive Discussion of the Possible Courses of Motion 
(Discussion Figures 27 and 28) 


2 

The same problem, although within a somewhat narrower scope 
(only system A), has already been treated by H. Bilharz in a report 
entitled "Ober eine gesteuerte eindimenslonale Bewegung" (On a 
controlled one-dimensional motion), Zeitschrift flir angewandte 
Mathe.matik und Mechanik, 22, pp. 206-215). The methods of investi- 
gation differ. H. Bilharz stressed particularly the mathematical 
representation, whereas we strove to emphasize the mechanical side 
of the problem. After the reports had been written independently, 
discussions held in common showed complete agreement of the results. 
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1. The Flip-Flop Control ^ (System Equation, Steering-Control Equation, 
Equation of Motion, Steering-Control Function, Systems of Q^ype A 
and B, Reversal Points of T^rpe a and b.) 

The controlled motions of a vehicle (which itself is regarded 
as & system of one degree of freedom.) are determined hy two equations: 
a) the equation of motion of the system ("system equation"), b) the 
steering-control equation. The system equation may be written, aB a 
rule, with greater or smaller approximations, in the 'form 


a'cp + b*q> + c’cp = Np 


( 1 . 1 ) 


where cp signifies the quantity characteristic for the motion 
(control or regulated quantity) — thus for the course control of a 
vehicle, for instance, the deviation from the course, p, the Bo-called 
adjustment quantity — for the course control the control— surface 
deflectionj the factor N is occasionally denoted as "control— surface 
effectiveness" . 


-^Recently a distinction in regulation technique has been made 
between the concepts "regulation" and "control" (Regulation Technique, 
Concepts and Symbols, edited by the 71)1 — Special Committee for 
Regulation Technique, outline, p. 3). A control is spoken of only 
when the variation in the position of the regulated quantity does not 
occur on the basis of a measurement of this quantity which is to be 
influenced; every variation of the regulated quantity based on a 
measurement of this quantity which is to be influenced is .called 
regulation. In this sense every so-called "automatic control of a 
vehicle" also is a regulation. In this case, however, the designation 
control is retained. In agreement with general colloquial usage in the 
technique of automatic vehicle control and particularly in aviation 
technique. For flip-flop control the outline mentioned above 
suggests the designation "open— shut regulator" or "in— out regulator" 
(elsewhere p. 83 ). Hie "system equation." (l.l) would be called, 
according to these suggestions, "equation of the- regulation distance", 
the control equation "equation of the regulator", etc. We limit 
ourselves here to these directives. 


If 
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If the regulation of the motion la performed by^ a control with,, 
continuous variation of the regulated quantity, the control equation 
reads in general 

k (k) n (V) 

21. v 3 =21 a y cp (1.2) 

K=0 V =0 


The number of the terms in the sums depends on the number of the 
sensed derivatives cp, q> ..... . and on the structure of the 
regulator (control machine) by which the number of the deriva- 
tives of p is determined. 

In (1.2) for instance the equation of the simple regulator 
position control 


P = + 


(1.2a) 


or of the simple flight— velocity control 


P = Sq9 + 


(1.2b) 


is contained (with use of two sensed quantities). 

By substitution of the equation (1.1) — solved with respect 
to p and differentiated as often as necessary — in the left side 
of (1.2), there results for the coordinate cp a hamogeneons linear 
differential equation of higher order with constant coefficients. It 
describes the controlled motion under Investigation. The roots of 
its secular equation for instance give information on the stability 
behavior. 

In this report a so-called flip-flop control (also called 
on-off -course control), in contrast to a control working continuously, 
will be Investigated. The characteristic of such a control consists 
in the fact that the adjustment quantity p does not vary continu- 
ously and proportionally to a linear combination of tha^ sensed 

quantities cp, <p, cp but can assume only two constant 

values +Pq or — Pq* The control equation reads in this case. If a 
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linear combination of the two quantities cp and cp, 
namely F = p^q) + Pg&j is sensed: 


P = — Pq sgn(p]Cp + P29) (1 . 3 a ) 

or 

0 = +Pq s S n ( Pi 1 ? + P2$) (1.3b) 


The function F = p^cp + Pg^j the sign of which, is in (1.3a) as well 

as in (1.3b) deciding the sign of the adjustment quantity g, will he 
called the (steering) control function below. 

If the control function F passes through zero, the adjustment 
quantity g is, in a discontinuous manner, "reversed" from, its 
positive value g Q to its negative value — g^ (or inversely) . Thus 

the points where the function F passes through zero are the "reversal 
points" of the control. As long as (like in this consideration) the 
zero passage of the control function 1 b decisive for the reversal, one 
may consider instead of 


F = p x «p + p 2 9 


the function 


F = q> + pep 


(1A) 


(where p = P 2 /Pi)* In this report the control function will always 
be used in the form (1.4), 

If p > 0, the control function F shows, with respect to the 
deflection, a lead by the time t T ; if p < 0, it shows a lag 

by t n (cf. fig. 1). 

Besides,. the time is counted so that from each reversal point 
a count starts with t = 0. The values cp (0) = cpQ and <p (0) = q>Q 
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are also called the reversal values. Because of F (0) = 0 the 
relation 


P fpo/^o 


(1.4a) 


exists between them. Let the length of the interval before the next 
reversal point be t^. 

A system, the control equation of which reads like (1.3a), is 
always denoted below as system Aj a system, the control, equation of 
which reads like ( 1 . 3 b), is denoted as system B. Thus one should note 
that the sign of the adjustment quantity 3 varies with the variation 
of the sign of the control function F, but that two possibilities of 
coordination exist; in system A one has • 


sgn 3 = -sgn F 


in system B, in contrast. 


sgn 3 = +sgn F 


Regardless of how the control-surface deflection may vary in the 
course of a motion, the correlation between the sign of 3 and- that 
of F is a priori determined by the apparatus. 

If one substitutes ( 1 . 3 a) or ( 1 . 3 b) in (1.1), the equation 


a*q) + b*cp + c*<p = +(sgn F)N 3 q 


or 


a’cp + b’cp + c’cp = +(sgn F)M _ ( 1 . 5 a) 


results as equation of the controlled motion of the system. In it 
the upper sign is valid for a system of the type A, the lower for one 
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of the type B. M = N&q is the value of the control moment * 


m = +(sgn F)M 


( 1 . 6 ) 


which stands on the right side of the equation (l.5&)j thus M is an 
essentially positive quantity. 

The variation of the control function F and of the control 
moment with time is shown in figure 2a for a system A, in figure 2h 
for a system B. Furthermore , a reversal point where the control 
moment m jumps from its negative to the positive value is called a 
reversal point a, one where the control moment jumps from the positive 
to the negative value, a reversal point h. According to the aforesaid 
it is clear that in the system A at a reversal point a the control 
function F passes through zero from positive to negative 
values (F < 0), at reversal point b, in contrast, from negative to 
positive values (F > o), whereas this correlation is reversed in the 
system B. 

In the form (1.5a) of the equation of motion four parameters, a*, 
b', c ' and M, appear ; furthermore all terms have the dimension of a 
moment (a generalized force). After division by the inertia factor 
(the moment of inertia) a f the equation (1.5a) assumes the form 

9 + 2Th9 + <u 2 cp = +(sgn F)b (1.5b) 


the number of parameters being reduced to .three, D, <n, b; now all 
terms of the equation have the dimension of an angular acceleration. 
The significance Is (as customary in oscillation theory) 


cu 2 = c’/a’j D = b’/2 \l a* c* and b = M/a* 


The upper sign on the right side again applies to systems of the 
type A, the lower one to systems of the type B. 

Occasionally instead of the form (1.5b) a further form of the 
equation of motion Is of advantage; it results from (1.5b) after 
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division by 


or 


and reads 


-%■ + 2D 2. + qj = +(sgn F) 4r (l-5c) 

of- CD £0^ 

Here all terms have the dimension of the coordinate 9 itself. 
Finally, with a new independent variable t = cot due to which the 
derivative becomes 


1 _d_ = JL_ 
£15 dt dT 


the differential equation (1.5c) could be changed to the form 

9” + 2Dp* + 9 = + (sgn F) \ (1.5d) 

of 

when the derivatives with respect to t are denoted by dashes. 

In all numerical calculations we shall further on put the 

parameter b/<x> 2 = M/e* equal to unity. This means that the scales 
used for plotting the results are at our disposal. The 

parameter b/o> 2 signifies the deflection (angle) 9=9^ that the 

control moment M produces on the restoring spring o'. This 
"static" deflection 9^ then Is the unit of the variable 9. Or, in 

other words, instead of the quantity 9, the new quantity 9 = 9/9^ 

is used as the dependent variable. Thus the differential equation 
ultimately obtains the final form: 

9" + 21)9* + 9 = +(sgn F) (I.50) 

Of the four parameters a*, b’, c’, and M appearing in (1.5a) 

or the three essential parameters D, cd } and b/05^ = M/e f appearing 

in (l-5c), two parameters, o> and b/cu , have been eliminated by the 
measures mentioned, namely, first, .different count of time (or. 
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respectively, use of derivatives «p/co, q>/oc , etc.)* second, having the 
scales of the variable cp and q?/co at disposal (or, respectively, 
different count of the deflection) so that as a sole parameter the 
damping coefficient. D remains. If nothing else is noted, the latter 
is put down for all numerical calculations as D = 0.1. 

For better maintenance of the illustrative quality, however, the 
equation of motion is not used in the form (1.5©); one stops at the 
form (1.5b) or (1.5c). 

Thus one calculates with the time t, using again as derivatives 

only quotients (p/co, tp/cu , etc.; one may also calculate in the 
coordinate cp itself, bearing in mind, however, that the latter is 
measured on a special scale, so that the unit denotes the static 
deflection cp^. 


2. The Finite Equation of Motion; the Equation for 
the Control Function and Its Derivative. 

The differential equation (1.5b) describes the motion in sections 
only. As soon - as the control function F changes itB sign, another 
differential equation appears, since the disturbance term in (1.5b) 
changes its sign. Likewise, the integral of the equation of motion, 
the time equation of the motion, applies for that reason only to one 
section between two zero— passages of the control function F; It 
reads 

Ant Xpt _ 

cp = A e + Be ^ + 


wherein 


= — 8 + IV = co — D 


= — 5 ■— IV = co — D 


(sgn F) -ky 

Oj 


( 2 . 1 ) 


+ i \Jl - D 2 J 

-i \/l - D 2 J 


( 2 . 2 ) 


■with 


8 = Dao and V = co Vl — D 2 . It has already been mentioned that 
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the count of the time in each interval is to start from zero. The 
time t = 0 denotes the beginning of an interval; a reversal point 
is located there. The derivative of (2.1) reads: 


cp = AlAe 


Xlt + X 2 Be X2t 


(2.3) 


Since for attainment of a real solution the integration constants A 
and B must he conjugate— complex, they may he written h_ 


A = Ce l€ , B = Ce -ie 


therewith one has obtained for (2.1) and (2.3) 


<p - 0 [e X l t+i£ + ? tign r) -i- 

Cl£ 


(2.1a) 


and 


<p = C 


A^t+ie Agt— i<Fj 

^2_ e + Age I 


(2.3a) 


A further useful form of the equation (2.3a) is obtained by 
expressing th e conjugat e— complex factors A^_ and Ag which* 

because of x/s^ + v 2 = m have the value <p, by this magnitude and the 
argument +or : 


Ap = aje icr 


A2 = cue 


(2.h) 


inversely, one then has 


= -*U_ Dj 

CD 


sin a 


-1- JT- 

o> v 


D 


COS <7 


(2.5) 
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Thua one obtains instead of ( 2 „ 3 a ) 


5> _ Q £ e A 'l t+i ( e+cr ) + e >^2'fc-i(e4 ; cf)| 


(2.3b) 


Analogously, the expression for the control function F = <p + p<p 
may be rewritten. Using (2.1a) and ( 2 . 3 b) one achieves, first 


t. rtT ^i^ + i e iff. lot— ie . — : icr~| 

F = C e x (1 + cepe ) + e r (l + cepe ) 


+ (sgn F) 

(XT 


( 2 . 6 ) 


If one expresses the conjugate— complex factors (in parentheses) by 
their magnitude C * and their argument ±t. 


1 + cepe icf = C*e iT j 1 + cepe icf = C*e 1 


(2.7) 


one obtains 


P 1.2/k+I ( e + ' r ) lot — i (e +t )~j _ . \b t r \ 

= CC*|^e x + e * 'J + (sgn F) -J- ( 2 . 6 a) 


therein the magnitude Is 


C* = \/l — 23*ep + o£p 2 


( 2 . 8 a) 


and the argument 


t = arctan 


1 - Deep 


( 2 . 8 b) 
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The derivative F of the control function shows the factors 
and Xg before the e— functions. They may again "be expressed 
according to (2.4) "by magnitude and argument. Thus one obtains 

= CC’[ © + e ^ ^ (2.^) 


Finally we indicate how the integration constants - C and e 
are related to the initial values cp Q and cp 0 /co: 



003 6 = ^[^0 ± ( B sn *) c) 

The equations (2.1a), (2.3b), (2.6a) and (2 .9) show the 

functions q>(t), 2&1, F(t), and i&l. We state them once 

more, splitting the coefficients X ± and X 2 into their real 
and imaginary parts: 
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cp(t) = Ce~^ t + e ~i(vt+e)J + ( B gn F) 

?(t) = Ce -5t j^i(vt+e+o) + e -l(vt+e+o)~| 


ar 


F(t) = CC'e^ [eKvt^) + e -i( v t4€+rjj + ( BgI1 F ) _L 

Elil a CC’e'"^' j^i(vt+e+ff+T ) + e — i(vt4s+cr+T )"j 
CD J 


The brackets may also be written in real form according to the 
de Moivre theorem: 


cp(t) = 2Ce cos(vt + e ) + (sgn F) 

CDr- 


= 2Ce~ cos(vt +e + c) 

03 

F(t) = 2CC , e — 1 ^ cos(/t + e + t) + (sgn F) 


CD 


= 2CC , e“' St cos (v t + e + a + t.) 


CD 


( 2 . 12 ) 


(2.13) 


In order to prepare for later applications, we Indicate explicitly 
what form the third equation of (2.13) assumes if the integration 
constants are expressed by the Initial values : 
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F(t) = + (sgn F) e -5t cos Vt 
cor 


i/l - D 2 


50 

to 


(1 - Hasp) + |flp 0 ± (sgn F) ~j(l>-top) 


e - ^ sin V t 
(2.14) 


+' (sgn F) — 

(XT 


At tlie reversal points t = 0 one obtains 


F(0) = 0 


and 


F(0) 

CJO 


9 0 /“ |_ 


& 


<p 0 

+ 2Dap,~ — + (sgn F) 
u o) " 


oV 


*o + < p 0 


_ i 




K 


(2.15) 


as one can recognize from F = tp + pep* using the differential 
equation (l.5t>). 


3. Bepresentation of the Course of the Motion. The Generating 

Motion. The Phase Plane. The Plane of the Eeversal Values. 

In order to represent the course of a motion one may use various 
expedients. The first, conventional one, consists of plotting the 
variable cp as a function of the time ; thus one obtains the 
deflection— time— diagram (q>-t— diagram) . The ordinate of this 
diagram, the deflection itself, may be obtained for special classes 
of motions (to which also belong those of interest here) as the 
projection of a "generating" motion. The generating motion of an 
harmonic oscillation is, for instance, a circular motion with constant 
angular velocity o>, a fact which is well known’ and often put to use . 
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The functions cp(t), Slil and. F(t) (cf. the equations (2.13)) 

CD 

all show the same damping factor; thus they may all he written in the 
form 


y = Ce -5 * 1 cos (v t + a) + y 0 (3-1) 


Such functions may, however, he represented as the projection of the 
motion of a point which moves with constant angular velocity v on 
a spiral (displaced from the zero point hy the distance y^) . Here 

the motion along the spiral is the "generating" motion. The spirals 
which lead to the functions cp, p/cr>, F, and f/cd are congruent. 

They differ only hy the angle at which the starting point cf the spiral 
lies (or, respect lire ly, hy which the spiral is rotated compared to an 
initial position) . 

The facts mentioned here are readily understandable. If one 
talces In addition to 


y = 


Ce 


— St 


cos(vt + a) 


(3-2a) 


the function 


x = Ce— 5 ^ sin(vt + a) 


(3-2h) 


and forms either 


z = y + ix or r = \/x 2 + y 2 

one arrives at the complex or real form of the equation of a 
logarithmic spiral In the form 


z 


Ce 


(— 5+iv )t+Ia 


(3.3a) 
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or 


r = Ce v ( 3 - 3^ ) 

Thus y(t) is one, x(t) the other rectangular projection of the 
motion of a point on the spiral. 

In the course of our calculations we have made wide use of this 
expedient of the generating motion on the spiral. The solution of 
transcendental equations which contain expressions of the form (3»l)j 
for instance F(t) = 0, is thus possible in a simple manner. 

The phase curve (curve in the phase plane) forms another means 
for following the course of a motion. Coordinates of the phase plane 
are the variables qp and cp or multiples of ^ them; for using variables 
of equal dimension we shall always select cp and <p/a>. The time t 
plays for such "phase curves" the role of a parameter; to each point 

of a phase curve there corresponds a constant value t. By differ- 
entiating the equation cp = <p(t) one obtains the second 
equation q>/o> = cp(t) /o>; both equations together form the representation 
by parameters of the phase curve. 

For a harmonic oscillation the equation of the deflection-time— 
diagram reads 


cp = A cos (cot + a) (3*^) 

the diagram represents a cosine— line. For the curve in the phase 
plane the parametric 'representation is , ^ • 

qp = A cos (cut + a) (3.4a) 

• 

— = -A sin(aot + a) (3.4b) 

a> 


Phase denotes in general physics any complex of the quantities 
determinant for a state (pressure, temperature...). A state of motion 
is completely determined by the coordinate and all its derivatives 

with respect to time; cp, <p', cp determine the phase of a state of 

motion (they determine a point in the phase space), cp and are two 
of the quantities determinant for the phase of .. a motion. They are the 
coordinates of the phase plane. In a problem. "of the second order" 
the phase is already completely determined by two quantities, cp 
and <p. 
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The equation of the curve, resulting by elimination of t, reads 

c p 2+ (I ) 2=a2 (3*^o ) 

it is a circle of the radius A. Thus to a harmonic oscillation a 
circle corresponds as generating motion as well as phase curve 
(fig- 3)- 

If a controlled motion of the kind named in sections 1 and 2 is 
described in the phase plane , all the reversal points for a certain 
system lie on a straight linej 
are made when F = cp + pq> = 0; 
count also starts anew* cp = cpQ 


or 

tan a 


the "line of reversals". Reversals 
however j since in that case the t ime 
and cp = <Pq. Thus one obtains 


p -_2a 

*0 


<P 0 

+ ~~~ = - cup 




(3-5) 


Thus the reversal points lie on the straight line, the slope tan a. 
of which has the value -afi. (if cup > 0, the straight lines go 
through the second and fourth quadrant, if cup < 0, through the first 
and third quadrant). 

The phase curves are similar to logarithmic spirals, but do not 
coincide exactly with them. Real logarithmic spirals as phase curves 
are obtained, however, if one places Into the phase plane Instead of 
a rectangular (rectilinear) coordinate system an oblique-angled one. 

This can easily be proved with the aid of figure 4. Let the 
angle between the negative abscissa— axis (-q>/tu — axis) and the 
positive ordinate axis (cp— axis) be the obtuse angle a as it Is 
determined as function of D by the equations (2.5). If one applies 
the cosine theorem to the shaded triangle, one obtains 


r 2 = 

cp + (sgn F) \ 

2+ (!) 2 - 8 

cp + (sgn F) ■— 

* 

or 


(XT 


(3-6a) 
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If cp and. cp/oj therein are expressed, as functions of _ the time 
according to the equations ( 2 . 13 ), there results 


— 2St 


P 2 

cos (V t + e ) + cos (vt + € + _cr) 


= 4c^e 

—2 cos(vt + e ) cos(vt + e + cr) cos p 


(3.6b) 


Splitting of cos (v t + e ) + a according to the addition theorem then 
yields L J 


r 


(2C sin cr) e 



thus the equation of a logarithmic spiral. 


(3-7) 


The lines of reversals j which have, as shown above/ in a rec- 
tangular coordinate system the slope tan a = -cup, have in the 
oblique-angled coordinate system mentioned the slope tan a = —tan t. 
This may be recognized with the aid of figure 4. For the angle a 
plotted there one finds 


tan a 


q> 0 sin a 
$q/cu — cpQ cos cr 


cup sin a 
1 + cup cos cr 


Hence, then results, with use of the equations (2.5) and (2.8) 

tan a = — = — tan t (3.8) 

1 - Hup / 

thus 

a = jr — t 


From equation (3*7) also follows that the angle subtended at 
the point of convergence of a spiral between the rays toward any two 
phase points has the value vt; thus it is a measure for the time 
required by the motion. Hence it follows directly that the angle 
between the rays toward the points of reversal has the value vt-j_ 

at the beginning and. at the end of the interval. 
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Since for system A the point of convergence of the spiral arc 
always lies on the side of the line of reversals turned away from 
the arc, for system B, however, on the side turned toward the arc, 
it follows immediately (cf . fig. 5) that for system A the angle 

is vt-j_ = «, for system B vtj_ = it (cf . also 5)* 

The fact that the phase curve in the oblique-angled system is a 
logarithmic spiral greatly facilitates following the course of the 
motion even through many intervals. One draws several windings of a 

- £ (vt) • 

spiral r = e v . Moreover one draws on a transparent paper 
provided with oblique— angled coordinate lines cp and cp/cu (cf . fig. 6) 
the points + b/cn 2 and the line of reversals. The line of reversals - 
for the points of which the control function is F(t) = 0 — then 
divides the phase plane into values pertaining to positive F(t) and 
into values pertaining to negative F(t). Let the starting point of 
the motion be P (cf. fig. 6). The control is to be handled 

according to system A. The control function F = cp + p<p has at 
the point P for the positive p a positive sign. ThuB the 

point -h/o^ must be fixed to the point of convergence of the spiral, 
and the spiral must be rotated so that it passes through P . The 

first point of intersection of the spiral with the line of reversals 
then is the initial reversal value Sq. At the reversal point S Q , 

F changes its sign. One must,- therefore, in continuing the 

construction of the phase curve, fix the point b/cu 2 to the point 
of convergence of the spiral; one may now draw the phase curve 
from Sq to the next reversal point S^. In this manner the entire 

course of the motion may be quickly constructed. 

In constructing such phase curves one observes easily two 
peculiarities which may occur in the course of the motion. First, 
in the proximity of the zero point of the phase plane there may 
exist reversal points starting from which the phase curve does not 
again meet the line of reversals (fig. 7)- Such a "missing" of the 
line of reversals may occur only in the case B, not in the case A, 
however strong the "damping" may be, for in the case A the point of 
convergence of the spiral lies "behind" the line of reversals. If 
the phase curve does not again meet the line of reversals, no 
further re7orsal takes place. Starting from the last reversal point 
the same functions cp(t) and (ip/u>)(t) continual l y describe the 
motion. ” The motion is a damped oscillation convergent toward the 
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point cp, = + -p-. This case will be treated more thoroughly; one 
a or 

speaks there of a "rest" of the control, the last reversal point being 
called "rest point". - 

The second peculiarity of the course of motion can also be 
easily understood from the phase pattern (cf. fig. 8). We have for 
instance a system A; we come from a reversal point S-,, the distance 

p 

of which from the zero point has the order of magnitude, b/co , to the 
farther Sg without having intersected the q>/a).— axis, and remain 

therefore at the left of the zero point on the line of . reversals . 

If we now continue with the construction in the customary manner, we 
find that the phase curve attempts to run further on the same side 
of the line of reversals, that is, the control function attained 
in S 2 the value 0, but did not, however, change its sign. Since, 

however, the "zero passage" is the condition for a reversal, the 
behavior of the control surface is not determined here. Only if for 
instance the reversal should occur somewhat after the zero passage 
of F (lagging of the control surface), we could construct further 
(cf. fig. 9)* One can see that the motion then assumes very high 
frequency and that its deflections decrease. The cases leading to 
such an undefined behavior of the "ideal" control are treated in 
detail in sections 4 and 7. 

In the further course of the investigation in addition to the 
phase plane another plane, showing a close relationship to it r will 
be used frequently. In the phase plane with the coordinates cp 
and <p/o) the motion may still be described continuously; a phase 
curve forms the sequence of all points of state. However, in our 
further investigations we shall often be content to consider instead 
of the succession of all points of state the se^uence ( of the reversal 
points, that is, of those values cp = cpQ and q>/co = cPq/cd which 

exist at the moment when the zero passage of the function F takes 

place. Thus from the points of the phase plane the reversal 

points cpQ, <Pq/cjd (which lie on a straight line, the line of reversals) 

are singled out and used for characterization of the motion. If one 
again forms a plane from the lines of reversals pertaining to all 
kinds of systems, one obtains the 9 q — cPq/o> — plane or_plane of the 

reversal values. In this plane of reversal values the motion of a. 
system is indicated by a sequence of points in isolated location. 

No further information is obtained concerning the course of the motion 
between the reversal points. The essential characteristics of the 
motion, however, — above all, the damping, the increase of the 
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amplitude , etc. — may be recognized even in such a representation, for 
the phase curves are logarithmic spirals with a point of convergence 

lying either at the point b/cu or — b/tu .. Thus one is certain that 
the motion between two reversal points does not show any peculiarities. 
The damping of the motion is assured when the sequence of the reversal 
points converges toward zero. 

However, whereas we placed an oblique-angled system of coordinates 
into the phase plane, we shall always place a rectangular one into the 
plane of reversals. The lines of reversals then again obtain the 
slopes tan a = -tup. This plane of the reversal points, 
the cpQ — 9 q/cjd — plane, will be discussed in great detail later on. 

An example is shown for illustration: Of the motion of a 

system represented by the phase curves of figure 10 (here drawn in 
Cartesian coordinates) only the reversal points appear in 
the cpQ — cPq/cd — plane (fig. 11). For the prescribed system with 

constant values p and cu the reversal points lie on the straight 
line of the slope tan a = -op . The points on straight lines of 
different slope pertain to other systems. However, it must be noted 
that the difference in the distances of two reversal points from the 
zero point is not yet a measure for the rate of damping. ThlB 
difference measures the damping per interval.. In order to arrive at • 
a damping velocity, the length of the intervals also must be taken 
into consideration (cf. section 6). 

A periodic motion is, for instance, represented by two reversal 
points repeated successively again and again (fig. 12) . However, 
these points need not by any means denote the maximum deflection 
which occurs during the motion; on the contrary, maximum deflections 
are reversal points only when i?q/cd = 0; such points lie on the 

ordinate axis. 


■4. Division of the Plane of the Reversal Points into "intercepts" 

and "Limit Points" ("Starting Points" and "End Points", 

Respectively), and into Reversal Points a and 

Reversal Points b. Respectively 

We start out considering the variation with time of the 
functions F(t) and m(t), as it is represented in figure 2. The 
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zeros of the control- function F(t) denote the reversal points. At 
these zero points the function F =9 + op ( 9 / 0 )) always shows a break, 

because cp/os^ ■ jumps (as follows from the differential equation, for 

instance ( 1 . 5 b)). The jump of 9 /a£ at a reversal point is 


— (+ 0 ) — - 2 - (- 0 ) = + ^ for reversal points a 


or 


or 


or 


(4.1a) 


■% (+ 0 ) — (- 0 ) = — for. reversal points b 

or- or or 

Because of the relation following from the definition of F 


t 

F 

Oi 


I 

= - + op 
0) 



there results, therefore, for the jump of F/co, which characterizes 
the break of the function F, at the reversal points a the value 


at the reversal points b 



F 

as 


(- 0 ) = - 2 op — p- 

or 


(4.1b) 


A positive value of the jump (right side of (4.1b)) for decreasing 
function F signifies a "breaking away from the perpendicular" 

(fig. 13a), for increasing function F a "breaking toward the 
perpendicular" if designations are used as they are customary in 
optics for the consideration of the ray path at the boundary of two 
media of different density. In contrast, a negative value of the jump 
signifies for decreasing function F a "breaking toward the perpen- 
dicular" (fig. 13b), for increasing function F a "breaking away 
from the -perpendicular" . "Breaking toward the perpendicular" increases 
the amount of the slope, "breaking away from the perpendicular" 
reduces it. 
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If one now considers what was said about the sign of the jump 
in the equations (4.1b),, together with what was said in section 1 
about the occurrence of reversal points a and' b in the systems A 
and B, the. following possibilities result: 


Table 4/1 



System A 

System B 

cup > 0 

breaking away from 

breaking toward the 


the perpendicular 

perpendicular 

cup < 0 

breaking toward the 

breaking away from 


perpendicular 

the perpendicular 


The correlations are independent of the "character" of the reversal 
points, whether a or b. 

As an example, we comment on the first case: If cup > 0, the 

jump values are, according to (4.1b), positive at reversal points a, 
negative at reversal points b. However, in the system A the 
function F decreases at reversal points a (where the jump value 
is positive) so that there a breaking away from the perpendicular 
occurs; at reversal points b (where the jump value is negative) it 
increases so that again a breaking away from the perpendicular occurs. 
The remaining three cases may be discussed correspondingly. 

From the fact that the break (change in slope) which the 
function F shows at the reversal points has a finite value, it 
follows that in those cases (that is, in quadrants where a breaking 
"away from the perpendicular" occurs) "directions of incidence" exist 
to which no "directions of refraction" on the other side of the time 
axis correspond; thus the function F does not again pass zero and, 
therefore, no more reversals occur. These cases correspond to the 
total reflections of the optical analogy (shaded angle section in 
fig. 13a). Inversely, in case of breaking "toward the perpendicular" 
"directions of refraction" exist to which no "directions of incidence" 
belong (fig. 13b). Thus in both cases no more reversal occurs; the 
zero points (reversal points) of the function F are no longer 
"intercepts" (of the controlled motion), but "limit points" where the 
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motion either must end ("end points") or may start ("starting points"). 
The directions of incidence without directions of refraction as well 
as the directions of refraction without directions of incidence lie 
in the angle sections shaded in figures 13a and h. The limitation 
of the angle- space is given by the condition that for the direction 
of incidence or refraction F = 0. If one substitutes for F the 
expression according to (2.15), one obtains 

(^f + 21 V ”0 * (sg^^o^o 2 - 0 (4 - 2) 

When D < 1, equation (4.2) represents for system A (upper sign) as 
well as for system B (lower sign) two ellipses in the plane of 
reversals, since sgn F may be either 1 or-. —1 (fig. 14) . If the 
third term has the negative sign, the equation describes the upper 
ellipse, if it has the positive sign, it describes the ‘lower ellipse. 
On the boundary of the ellipses lie the "limiting cases of total 
reflection"; within the ellipse lie limit points; outside of it, 
intercepts. - 

Whether the limit points lying in the interior of the ellipse 
are end points or starting points may be decided for instance in the 
following manner: For end points, breaking away from the perpen- 

dicular always occurs; for starting points, breaking toward the 
perpendicular. The table 4/1 , which classifies the quadrants in 
this respect, gives, therefore, information also about the question 
where end points and where starting points lie in the space Inside 
of the ellipses. The results of this consideration are plotted in 
figures l4a and l4b. The ellipses represented by equation (4.2) are 
occasionally called the "small limiting ellipses" to distinguish 
them from other curves to be discussed later. 

As explained in section 1, the reversal points a and b 
differ by the sign which the function F shows at the point zero. 

The function F(0 )/od is given in (2.15). Thus F(0 )/cd is positive 
as long as both factors, and the brackets K, show, the same 

sign; f/cd is negative when the factors have 'different "signs . The 
limiting case f/cd = 0 occurs when the brackets vanish; the 
significance of this case was discussed Just now. The equation K = 0 
represents the "small limiting ellipses" discussed above. Eeversal 
values cPq, 9q/o> from the space inside of the ellipses lire 



HACA TM 1237 


25 


for K < 0; reversal values from the space outBide of the ellipses 
are for K > 0. ‘ 

We inquire first: What regions of the cp Q — cp 0 /os — plane , for 

the system A, contain reversal points a. To system A pertains the 
upper sign of the third term in K (2.15) • Furthermore, one has fear 
every reversal point a (cf . fig. 2a) in the system A (after the 
reversal) sgn F = — 1. Therewith the term in question assumes the 
negative, sign; thus of the two "small limiting ellipses" only the 
upper one is considered.' Reversal points a in the system A require, 
moreover, F-^-0, thus different signs for the factors and K. 

In the right semiplane where cPq/cd > 0, K must he less than 0; the 

reversal points a lie, therefore, in the space inside of the (upper) 
ellipse . In the left semiplane where (Pq/o 5 < 0, K must he greater 

than 0; the reversal points a lie in the space outside of the 
same ellipse (cf. fig. 15a). 

If one investigates reversal points h, one obtains sgn F = 1, 
the third term in K becomes positive, the lower one of the limiting 
ellipses is of importance. Since F /ca must now he greater than 0, 
the reversal points h lie in the right semiplane in the space 
outside, in the left semiplane in the space inside of the (lower) 
ellipse (fig. 15b). 

The investigation of system B is similarly found. Here the third 
term of K (2.15) has the lower sign. For reversal points a, 
sign F = 1, so that the negative sign for the term remains; by this 
fact (as for reversal points a of system A) the upper ellipse applies. 
The reversal points a require that F/cu > 0; they lie in the right 
semiplane in the space outside and in the left semiplane in the space 
inside the (upper) ellipse. (Fig. 15c.) For reversal points h, 
because^ the sign F = —1, the lower ellipse is of importance; 
since F/to <0, the reversal points b lie in the right semlplane in the 
space inside and in the left semiplane in the space outside the (lower) 
ellipse. (Fig. 15d.) 

If one now visualizes the figures 15a and 15b of system A 
superimposed, one can see that within the ellipses at the upper left and 
lower right no reversal. points, either a or b, are situated. 
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These points we re.cognized In section 3 as end. points of the motion. 
In contrast, the starting points which lie within the ellipses at the 
upper right and lower left are included in the consideration, since, 
the definition of the reversal points makes use of the sign of the 
function F after passing zero, and this sign actually is explained 
for starting points. It even becomes evident that these last-named 
regions are doubly covered over when figures 15 a and 15 b are 
superimposed, and thus contain points of type a as well as of 
type b. The motion may, therefore, begin from all. starting points 
In two different manners. Outside of the- ellipses lie intercepts 
which afe now unequivocally divided into reversal points a and b. 


5 . Curves of Constant Interval-Length Vt-^ in 
the Plane of the Reversal Values (isochrones) 
and Their Particular Cases 

The length of the intervals Vt is different, for the separate 
reversal points. We Investigate here the geometrical loci of all 
the reversal points for which the following interval is of equal 
length. It can be shown that these geometrical loci are ellipses. 

From the condition that 


F(t-L) = 0 


(5.1) 


follows with the aid of equation (2.14a) after a few transformations 
which may easily be supplemented 


(%f ± T) 4r 


— Vt. • 

I 5 cos vt. — e v x 

D + 1/1 - I r — — — l 

sin vti 


2a 

CD 


(5.2) 


+ 2Dcp o % ± ^ F) \ <P 0 + ?o 2 = 0 

CD (Xf- 


For D < 1 this is the equation of an ellipse in the plane of 
reversals with-.Vt^ as a parameter. 
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Tiie numerical investigation shows that for the system A only- 

ellipses may appear for which Vt-j_ = Jt, for the system B only 

ellipses for which vt^ = jt. This fact is immediately understandable 

if one considers that in the system A the control moment essentially 
strengthens and in the system B essentially weakens the restoring 
moment. A glance at the phase curves of section 3 (fig. 5) also 
shows that the length of the interval compared to an uncontrolled 
motion is diminished in the case A, increased in the case B. 

In section 3 it had also "been recognized that, if the length 
of the interval t^ is increased, the phase curves finally no 

longer intersect the line of reversals so that no further reversal 
occurs. The last reversal point is called "rest point". In the 
deflection— time-diagram the limiting case of such a motion appears 
as indicated in figure 16. This limiting case is characterized by 
the additional condition (besides ( 5 * 1 )) 


F(t3_) = 0 


(5-3) 


This equation (5-3) singles out solely one ellipse from the 
family of ellipses and determines the parameter t-j_. 

Using equation (2.14) and its derivative, one obtains from the 
conditions ( 5 . 1 ) and ( 5 . 3 ) the transcendental relation 


5 -®(vbi) 

cos vt^ — — sin vt — e =0 


( 5.*0 


for the parameter vt^. It is satisfied 


1) for Vt = 0 


2 ) for a value different from zero vt., = vt n 
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This last value is the one that corresponds to the "limiting rest 
condition". _ Thus an ellipse is determined in the plane of reversals 
by vt-j_ = vtj_j it is occasionally called the "large limiting ellipse". . 

The maximum interval length following from equation (5*4) has., 
for D = 0.1, the value 


vt-j_ RJ 3OO 0 


Such rest points can occur only in the system B. Besides, it can 
easily he seen that the "small limiting ellipses" mentioned in 
section 4 also belong to the family of ellipses of constant interval- 
length. Their parameter has the value v t^ = 0; one recognizes at 

once that one obtains from equation (5-2), by performing there the 
limiting process vt-^ ->0, the equation (4.2). 

As mentioned at the beginning, the amount of damping is given 
the value D = 0.1 in all numerical calculations. All figures are 
drawn accordingly. The figures 17 a and 17b show (for reversal 
points a) how the small and large limiting ellipses appear when D 
assumes other values. The curves are plotted for D = 0.1 
and D = 0.5* In any case the curves go through the points 0 and 1. 
For D = 0 the small, and likewise the large, limiting ellipse turns 
into the circle through the points 0 and 1 symmetrical to the 

If one places — as we did — into the plane of the reversal values 
a Cartesian coordinate system, the curves of constant interval- 
length t^ are ellipses. Without explicitly carrying out the 

transformation we want to note that all these ellipses become circles 
if one applies instead of the rectangular the . oblique— angled 
coordinate system that was used in section 3 for the phase planes- 
(where the negative q/a> —axis forms with the positive cp —axis 
the angle a). The equation (5*2) then is the equation .of the 
family of these circles; the angle a appears if one introduces it 
according to (2.5) instead of the 'amount of damping D. . . 



6. Energy Consideration and Periodic Solutions 

For the present problem the flip-flop control has the purpose 
of achieving a much stronger damping of the motion than attained by 
means of the existing damping force alone. However, in such a 
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controlled system not only damping, tut also periodic and amplitude- 
increasing motions are possible since the control moment may not 
only remove energy from the system hut may also supply energy to it. 

In order to obtain a survey of where in the cp^ — — plane 

lie the reversal values for damping and where for periodic and 
amplitude— increasing motions, we inquire first about those points 
In the plane of reversal values from which motions start that are 
deprived of energy or supplied with it by the control moment. 

The energy supplied to the system by the control moment m Is 
given by the expression 

^1 

U z = /* mp dt (6.1) 

0 

With 

m = + (sgn F) M (6.2) 


(the upper sign again being valid for systems A, the lower for 
systems B) It becomes 


U z = + (sgn F) M (cp 1 - cp Q ) (6.3a) 

Ve now investigate what sign this expression has for reversal 
points in the individual quadrants. We consider first system A. 

In section h we determined that for system A reversal points a lie 
in the second and third and reversal points b in the first and 
fourth quadrants. With reversal points a in system A sgn F = — 1; 
with reversal points b sgn F = 1. Under the further assumption 
that the course of the motion is always such that cpQ and cp-^ lie 

on different sides of the zero point^, the factor (cp^ — cp^) becomes 


5 

"Tor intervals leading to end points or starting from starting 
points, cPq and <p^ may happen to have the same sign. In that case 

one must consider them separately; however, the result will again be 
the one shown in table 6/1. 
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negative for motions starting from reversal points in the first and 
second quadrant, positive for motions starting from reversal points 
in the third and fourth quadrant. " 

Talcing these facts into consideration one finds the first column 
of the following table. 


Table 6/1 


sgn U z 


Motion starts from a 
reversal point in 

System A 

System B 

1 . quadrant 

1 

-1 

2 . " 

-1 

1 

3- 

1 

-1 

k. 

-1 

1 


For the system B first the first sign in ( 6 . 3 a) is reversed. 

Furthermore the correlations of (sgn F) to the reversal points a 
and b are interchanged; however, since the reversal points a 
and b simultaneously vary their position in the quadrants, the 
correlation of (sgn F) to the quadrants is maintained. Likewise the 
correlation of sgn (cp^ — cp^) to the quadrants remains intact. That 

means the sign of U_ in system B is, for all quadrants, the opposite 

of that in system A (column 2 of the table 6/1). 

Where sgn U z = 1, energy is supplied to the system by the 
control moment, where sgn U z = — 1, energy is removed. However, 

where energy is removed by the control moment . as. well as. by the 
damping moment, the controlled motion certainly is damped even more 
strongly than the uncontrolled one. Where, on the other hand, energy 
is added, it is damped less and the amplitude may even become constant 
or Increase. 

The motions starting from reversal points, in those quadrants 
where the control moment supplies energy must now be considered more 
thoroughly in order to determine the position of the reversal points 
from whence start periodic, damped, or amplitude— increasing motions. 
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The energy supplied, "by the control moment is given "by 
equation (6.3a); it may "be put in the form 


U z = c’ 


+ (sgn F) 


J? (<p ! - <p 0 ) ] 


For the energy consumed, "by the damping force one has 


U v = 


(b'tp) 9 dt = a’2S ^ cp 2 dt 


(6.31) 


(6.4a) 


If one introduces 9 according to equation (2.3b) and integratesj 
one obtains 


U y = c ’DC 1 


: j^ e 2X l t l + 


(2€+cr)i 2X2t]_ — (2s+o)i 


+ e 


(2e+a)i — (&+c)i 0 ~25t« o 

— e — e — — e ^ — 

D D_ 


(6.4b) 


If one expresses the e-functions "backwards again "by 9^ 9 q/<x>j 9-^ 
and 9^/oij one obtains 


U_ 



+ 2C 2 (l-l£)(l - e^ 8tl ) 


T (sgn I) d(S£ - ®L) 


( 6 . 5 ) 
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and, accordingly , for the difference of the two energy values 


& = u v - u z = 


1 1 
= ("b f q> )cp dt — J' nrp dt 


c, «»i^-’’ov) + 2c2(1 


j + 2C 2 (1 - T?) (l - s' 26 * 1 ) 


~2Sti\ (6.6) 


+ (sgn F) i £>(%■ - h) + („ 0 _ ^ 


First, it shall he shown that vanishes for reversal points 

belonging to periodic motions. 

The periodic motions themselves may be found without an energy 
consideration: For periodic motions there exist between the reversal 

values cpQ, <Pq/o> at the beginning of the interval and the reversal 

values cp^j q>j_/cD at the end of it the relations 


and 


90 = -c Pl 

% = _ < h 

CD • CD 


(6.7a) 


(6.7b) 


The period then is the double interval— length T = 2 t^ . With the 

aid of the equations (2.13a) and (2,13b), (6.7a) and (6.7b) may be 
put in the form 


— Sts 


2 C cos e + (sgn F) — = - 2 C cos(e+vt 1 ) e + (sgn F) — (6.8a) 

co d • CD 2 * ; 


and 


— &ti 


cos(e + a) = - coa(e + a + vt^) e 1 


(6.8b) 
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The equations ( 6 . 8 a) and. ( 6 . 8 b) represent two relations between the 
integration constants C and e and the quantity a characteristic 
of the system on one hand and the interval length vt^_ on the other. 

If one expresses the angles e and a with the aid of (2.11) and 
( 2 . 5 ) by the reversal values cp^ and ^q/cd, one achieves (after a 

longer calculation) the equations 


(\A - d 2 ) 8111 ^ 1 - 3)2 - 


Bin Itct., 


cp Q = + (sgn F) ~ 

40 cos(\/l — D“ cot^J + cos Eustj 


^ = ± (sgn F) 4 

CD (vf 


sin ( \/l - D 2 CDt-J 


\fT~& cos^l — D 2 cot-^ + cos Itot^ 


( 6 . 9 a) 


( 6 . 9 b) 


as parametric representation of the "curve of the periodic solutions" 
in the plane of reversal values. 

This curve is a spiral— like curve (cf. figures 28a and 28b) ; the 
interval length t^_ is a parameter. The curve starts for t-j_ = 0 

at the origin of the coordinate system and runs for t^— toward 
the point <Pq = ±1, ?q/cd = 0. For Vt^ = it the curve intersects the 
ordinate axis. The points vt^ < it belong (according to the 
aforesaid) to the system Aj the points vt-j. > it to the system B. 
Moreover , the curve ends where v t-^ becomes vt-^. It ends* therefore^ 

on the "large limiting ellipse" (which limits the region of the rest 
points) ; it can be shown that that point is simultaneously the point 
of intersection with the small, li mi ting ellipse. 

If one now substitutes in (6.6) — after having expressed C 
also by cp^ and <Pq/cd — the reversal values according to ( 6 . 9 ) and 

takes into consideration that moreover cp^ = -cp Q and cp^/rn = -cp^/o^ 

one finds 


or 


£U = 0 
U v = U z 


( 6 . 10 ) 
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Aa before, the value D = 0.1 of the amount of damping is taken 
as a basis also for the curves of figures 27 and 28. For other values 
of the amount of damping the curves of the periodic motions appear as 
indicated in figure 18 which refers to reversal points a. For D = 0 
the curve degenerates into a straight line which coincides with 
the cPq/o) — axis. 

Furthermore, the sign assumed by the energy difference 
outside and inside of the curve of the periodic motions (6.9) is 
of interest. “ 


A generally valid discussion of the sign with the aid of the 
equation (6.6) seems difficult. However, a large number of numerical 
examples show that for system A in first and third quadrants, for 
system 33 in the sector between the cpQ — axis and the ray OSQ (fig. 28) 

the motions always converge toward the periodic ones, from reversal 
points outside of the curve of periodic motions with decreasing 
deflections and from reversal points inside of that curve with 
increasing deflections. In the sector between the rays OSQ and OR 
the motions go from "inside" points with increasing deflections toward 
the periodic ones. 

The arc RS of the curve of the periodic motions is "unstable", 
that is, from reversal points in its vicinity the motion goes, in 
case the point lies "inside", toward a periodic motion on the arc 
section lying farther to the outside, in case the point lies ’’outside", 
the motion is damped in the customary manner (leading to a rest point). 


A good survey of the course of the motion discussed just now and 
therewith of the problem of the removal or addition of energy may be 
obtained from figures 19 to 2b. Figures 20, 21, and 2b show, against 
the non-dimensional interval length V tp the pertinent value of 


reversal cp^ at the beginning and 


at the end of the interval 


plotted for various values p. If one selects a first reversal 

point, for instance in figure 21, all the following are easily determined 

by drawing in a stepped line. Starting from the reversal point S^, 


one no longer obtains an intersection point with the curve cp Q . 


The 


point corresponds to a point in the domain of the small limiting 


ellipse. For p = 0 one must plot 9 q/o> since all cp n = 0, 
(n = 1, 2, 3, ...), see figure 19. For p-> + »,vt^ always 
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equals it and, moreover., <Pq always equals zero. Thus one has to 
select here also another manner of plotting. We selected the one 

used hy K. BBgel^ since for p— our problem transforms Into the 
one treated by B&gel. For this case a simple relation between the 
reversal values cp Q and cp^ may be given (cf. figs. 22 and 23). 

A further problem arises in connection with the energy consider- 
ations. Even the free oscillation of a damped system causes removal 
of energy. If the control of the system is to damp the motion more 
rapidly, the decrease of the maximum deflections must be stronger 
than for the free oscillation. For the free oscillation the ratio 

— — it 

of two successive maximum deflections is e v . For the 

controlled oscillation the ratio is not constant, and the time t w 

M 

between two maximum deflections is not constant, either. It is, 
however, assumed that t^ always lies between the correlated 

interval t 1 and’ t.. (cf. fig. 25a). This relation may be proved 
L I ± I1 

with the aid of the phase diagram (fig. 25b). For this purpose the 
construction is not continued from the reversal point Sg in the 

customary manner, but S2 is — since we are interested only in the 

absolute values of the maximum 9 — values — mirrored at the zero 
point; then one may continue to operate in the second interval with 
the same point of convergence, as can .be understood immediately. One 
has ^SqPS! = Vt-j^j <5S 1 , PS 2 ‘ = Vt^; = vt M . By 

construction. 


< $ S-lPS^ < SQpSi 

By the construction a simple correlation between t^ and the 
adjoining interval lengths is given. For this reason there iB for 
system A vt^ ^ it and for system B it = Vt^ = Vt z (cf . p. 27 ). 


^B&gel, E. : Das Yerhalten ged&mpf ter und auf schaukelnder freier 
Schwinger unter der gleichzeitigen Einwirkung einer konstanten 
Reibungskraft . (The behavior of damped and amplitude— increasing free 
oscillators under simultaneous influence of a constant frictional 
force.) Ing—Archiv 12, p. 2^7. 




36 


NACA TM 1237 


One may now plot the ratio of two successive maximum values 
(construction cf. fig. 25b) against vt^ (fig. 26). The_ uncontrolled 

oscillation is represented in this diagram hy the point k (vt^ = )t). 

In order to he enabled to compare the damping of the motion in the 
controlled and in the uncontrolled system, one has to take into 
consideration that in the controlled system the change of the maximum 
deflection occurs at different times. The rate of damping in the 
controlled system is (* of the uncontrolled motion was denoted 

hy of the controlled motion hy 


_ _ 9Gn+ 1 ~ ^Gn 
G ~ Vtw 

In the uncontrolled system, that is, for a free damped oscillation 

IvUn+ll < P[Jh I 
""a 

Only when the ratio v q/ v jj exceeds 1 does the control fulfill its 
purpose. Thus one forms the ratio y q/ v jj for e/jual Initial 
value q^ = q^: 


L 

r a 

VGn+1 

jvGnl 

l-i 1 9 Gn.+l | 

I^Gnl 

Y u “I 

%n+l 

/<!>Uh/ 

1- i vtM ( 

( r \ Vt M 
\® v 1 ~ 1)2 ~ l ) 
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The numerator 


The denominator \e 


lw / 

i * -** 


I^Qil 


can immediately "be taken from figure 26 ^ . 


D 

Vl - D 2 


- 1 


Vt 


M 


TT 


represents a straight line 


which is drawn in for comparison. It goes through the zero point and 
through the point JE. For all points lying "below the drawn— in 
straight line f , the ratio is y q/ y jj > that is, the controlled 

motion is damped more quickly than the uncontrolled one. This, 
however, signifies that, as a rule, the system A will "be preferable. 
We shall return to this problem In section 7. 


7- Comprehensive Discussion of the Possible Courses 
of Motion (Discussion of Figures 27 and 28) 

The perceptions gained in the previous paragraphs are summarized 
in figures 27 and 28. In both figures the type of line used indicates 
the type of the reversal points: Solid curves correspond to the 

reversal points a, dashed ones to the reversal points b. Figures 27a 
and 27b show for system A and for system B, respectively, the curves 
of constant interval length tp (isochrones) . Here one recognizes 
the fact mentioned In section 5, that in the system A the parameter 

is t^ =" jt/v , In the system B, in contrast, t^ "S’ n/v . To the family 


^Regarding the limit rest In figure 26 the following remark must 
be made: On p. 35 eld- upper limit is determined for the time 

interval vtj^ between the successive maximum values. Only maximum 

values within the sphere of normal operation of the control are 
considered, for Instance cp^ and cpQg in the sketch, figure 26, 

upper right. If, however, in the system B rest occurs in the course 
of a motion, one is inclined to select (in considering the damping of 
the controlled system) as the last pair of maximum values 

and The time interval ( v between these maximum values, 

however, may be larger than Vtp, as can be easily taken from the 
phase curve: « < ( v "* : 'm)r < 
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of these curves of constant interval length also belong those 
for t-^ = 0 and t^ = t^, which had been denoted as the 11 small" 

and "large" limiting ellipses in sections 4 and 5, and the properties 
of which had been discussed there. In figure 29 the small and the 
large limiting ellipse are (for the system B and reversal points a) 
drawn separately once more and the significance of the individual 
domains resulting from section 4 is indicated. 

Besides the small limiting ellipses which enclose the domain of 
the limit points, those domains are marked in the two figures 28a 
and 28b, whence one arrives at an end point after one or after two 
intervals. The domain whence a single interval leadB to an end point 
is subdivided once more and it is noted whether this end point liee on 
this or the other side of the zero point. 

For system B, aside from end points, rest points of the motion 
are also possible. Reversal points in the crescent— shaped domains 
lead after one or after two steps to such rest points. 

After all, that has been said so far, the course of the motion 

adjoining any reversal point is now clear and thus the time has come 

to decide on what parameters are useful for a desired course for given 

initial values. In section 1 it was shewn that, when the actual 

time t is used as independent variable, the equation of motion 

contains the parameters D, oa, b, and p. Let us assume that D 

and cs are prescribed, that is, that we are__dealing with an 

oscillating system the damping and natural frequency of which are 

given. Let the damping be insufficient and the requirement be made 

that the damping be improved by application of an intermittent 

control. It is, therefore, desired to attain in the shortest possible 

time, for prescribed initial values 9 and 9 , negligibly small 

8 . 0 . 

values 9 and 9. We may now select b/cu and p and decide, 

moreover, on control system A or B. If, for instance, a large b/cu^ 
is selected, it may happen that almost all reversal values pertaining 
to our maximum possible initial values fall into the domain of the 
small limiting ellipse. That would mean that on the basis of the results 
of this report ("ideal" flip-flop control without lag) the motion becomes 
undefined after a short time or increases toward a periodic motion 

(cf. fig. 14 and p. 3^) • If,, on. the other hand, a small b/co 2 is 
selected, the reversal values in question will lie within a domain of 
their plane which is sq large that the limiting ellipses and the two. 
curves of the periodic motions lie entirely within it. If now for 
instance a negative p is selected, we would have in system A for a 
large domain on initial reversal values an increase of amplitude; for 
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positive p, on the other hand, in general, damping extension into 
the domain of the small limiting ellipse would occur. Moreover it 
would have to he investigated whether the rate of damping is 
sufficiently high. If system B is selected, there would occur, for 
negative p, decrease down to the rest ellipse; for positive p an 
increase of amplitude occurs for a large domain (within the curve of 
the periodic solutions). Since the ellipse of the rest points is 
always larger than that of the end points (small limiting ellipse), 
system B with negative p will certainly he more unfavorable than A 
with positive p. Let us, therefore, decide on system A and 
positive p. Now arises the problem what magnitude of p to select. 
For this purpose one must consider the curve of damping against the 
time t, that is, the figures 19 to 22 and figure 26. According to 
these figures, one is at first very much inclined to select a rather 
large p. However, the following fact is opposed to that: As we 

know, one enters rather soon the domain of the mnan limi ting ellipse 
in which the motion, under the effect of a control without lag, is no 
longer defined. It has been pointed out before, on p.20, that one 
again obtains defined phenomena of motion in this domain, if one 

Q 

takes the lag into consideration (cf . the report by K. Scholz 0 made 
in connection with our report) . It is now shown that, with these 
facts taken into consideration, smaller values of p are preferable 
since they cause a faster damping of the motion within the small 
q 

limiting ellipse . 


Translated by Mary H. Mahler 
National Advisory Committee 
for Aeronautics. 


o 

°Scholz, K, : tJber die Bewegungen eines Systems von einem Frei— 

heitsgrad unter dem Einf luss einer Schwarz-Weiss— Steuerung mit 
Schaltverschiebungen. (On the motions of a system of one degree of 
freedom under the effect of a flip-flop control with lag of 
reversals . ) 

^Frof . Fischel (BFS) called our attention to the fact that it is 
shown in a report by Gelling (shortly to be published) that ior very 
small time lags the motion turns into a creeping toward 

zero (qp = e — fc /Pj. 
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(a) 


(b) 



Figure 2.- Course of the control function F and of the control_mornent m 

for system A and system B. 



Figure 3.- Phase curve of a harmonic oscillation. 
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(a) (b) <c) (d) 

Figure 15.- Domains of the reversal points a and b. 
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Figure 19.- Survey of the decrease and increase, respectively, of the reversal 

values of controlled motions. 



Figure 20.- Survey of the decrease and increase, respectively, of the reversal 

values of controlled motions. 



Figure 21.- Survey of the decrease and increase, respectively, of the reversal 

values of controlled motions. 























Figure 27.- Curves of constant Interval length h (Isochrones) in the plane of the reversal values 

CrrcfoTn A » — 
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Figure 28(a).- Division of the plane of reversal values. System A 
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System B 
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Figure 29.- Large and small limiting ellipses. System B. Reversal points a. 
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